Quantum renormalization group approach to geometric phases in spin chains 
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A relation between geometric phases and criticality of spin chains are studied by using the quantum 
renormalization-group approach. We have shown how the geometric phase evolve as the size of 
the system becomes large, i.e., the finite size scaling is obtained. The renormalization scheme 
demonstrates how the first derivative of the geometric phase with respect to the field strength 
diverges at the critical point and maximum value of the first derivative and its position scales with 
an exponent of the system size. This exponent is directly associated with the critical properties 
of the model where, the exponent governing the divergence of the correlation length close to the 
quantum critical point. 
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I. INTRODUCTION 

Quantum phase transition (QPT) has been one of the 
most interesting topics in the area of strongly corre- 
lated systems [l|. at zero temperature the properties 
of the ground state may be changed drastically shown 
as a non-analytic behaviour of a physical quantity by 
reaching the quantum critical point. This can be done 
by tuning a parameter in the Hamiltonian, for instance, 
the magnetic field or the amount of disorder. Tradi- 
tionally such a problem is addressed by resorting to no- 
tions like order-parameter and symmetry breaking i.e., 
the Landau-Ginzburg paradigm [2]. In the last few years 
a big effort has been devoted to the analysis of QPTs from 
the perspective of Quantum Information the main 
tool being the study of different entanglement measures 
[rOi] . In view of some difficulties [ll| , attention has shifted 
to include other, potentially related, means of character- 
izing QPTs fl^l. One such approach centers around the 
notion of geometric phase (GP). GP have been offered as 
a typical mechanism for a quantum system to keep the 
memory of its evolution in Hilbert space. Such phases 
were introduced in quantum mechanics by Berry in 1984 
[Tsj . Since then, geometric phases became objects of the- 
oretical and experimental research [l4j, whence it was 
noticed that they are related to a number of important 
phenomena in physics (Tsl, s uch as Aharonov-Bohm fld\ 
and quantum Hall effects [l7[. In recent years, this inter- 
est has increased due to their applicability in quantum- 
information processing [l8|. In other words, GP becomes 
increasable for product states of composite systems since 
the uncorrelated subsystems pick up independent geo- 
metric phase factors. However, GP are induced by quan- 
tum entanglement, if the full state is pure. On the other 
hand, classical correlations and quantum entanglement 
can coexist in mixed quantum states, which means that 
the forms of mixed state geometric phases in [l^, ap- 



plied to the path of relative states, may contain portions 
from both these types of correlations. Nevertheless, their 
connection to the quantum phase transitions has been 
manifested recently in Ref . |20l - [23 ] . where they has been 
shown that the geometric phase can be used to investi- 
gate the critical properties of the spin chains The 
critical exponents can be evaluated from the scaling be- 
havior of geometric phases [2l| . Therefor the geometric 
phase can be considered as a topological test to manifest 
quantum phase transitions ^2^. These general relation 
roots at the topological property of the geometric phase, 
which describes the curvature of the Hilbert space, and 
especially has direct relation to the property of the de- 
generacy in quantum systems. The degeneracy in the 
many-body systems is critical in our understanding of 
the quantum phase transition. Thus the geometric phase 
is another powerful tool for detecting the QPT. 

Our main purpose in this work is to use the ideas of 
quantum renormalization group (QRG) [23] to study the 
evolution of the geometric phase of the spin models. To 
have a concrete discussion, the one dimensional S = ^ 
Ising model in transverse field (ITF) has been considered 
by implementing the quantum renormalization group ap- 
proach 0-11, milBl- To the best of our knowledge, the 
GP properties of the exactly solvable models have been 
studied up to now. This is the first work which addresses 
how to get GP properties of the models which are not 
exactly solvable using QRG. However we show that in- 
vestigating the critical properties of models by studying 
the GP properties of models using the QRG method is 
more convenience and accurate than its application on 
entanglement (concurrence) . 



II. THEORETICAL MODEL 

We have considered the ITF model on a periodic chain 
of N sites with Hamiltonian 
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where J > is the exchange couphng and A is the trans- 
verse field. From the exact solution [2^ [23| it is known 
that a second order phase transition occurs for Ac = 1 
where the behavior of the order parameter or magneti- 
zation is given by < >= (1 — A)^/^ for A < 1 and 
< cr^ >= for A > 1. 

III. QUANTUM RENORMALIZATION GROUP 

The main idea of the RG method is the mode elimina- 
tion or thinning of the degrees of freedom followed by an 
iteration which reduces the number of variables step by 
step until reaching a fixed point. In Kadanoff's approach, 
the lattice is divided into blocks. Each block is treated 
independently to build the projection operator onto the 
lower energy subspace. The projection of the inter-block 
interaction is mapped to an effective Hamiltonian [H^^^) 
which acts on the renormalized subspace [23 - l26| . 

To implement QRG the Hamiltonian is divided to two- 
site blocks 

N/2 

1=1 

and the remaining part of the Hamiltonian is included in 
the inter-block part 

N/2 

1=1 

where (t"j refers to the a-component of the Pauli matrix 
at site j of the block labeled by /. The Hamiltonian of 
each block {hf) is diagonalized exactly and the projec- 
tion operator 

Po = |V'o)(V'o| + |^i)(V'i|, 

is constructed from the two lowest eigenstates. Where 
Itpo) is the ground state and \tpi) is the first excited state. 
In this respect the effective Hamiltonian 

is matched to the original one (Eq.(IT])) replacing the 
couplings with the following renormalized coupling con- 
stants. 

J' = J-^, q = X+ VA2 + 1, A' = A^. (2) 

1 _|- /-/-^ 



IV. GEOMETRIC PHASE AND 
RENORMALIZATION GROUP APPLICATION 
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FIG. 1: (Color online) Evolution of the geometric phase under 
RG versus A. 



be described by applying a rotation of around the z 
direction to each spin |2Q|, i.e., 

N 

il4> = 90 = n expl-ifZ-crf /2). 

The critical behavior is independent of as the spectrum 
of the system is independent [21*1 . The geometric phase 
of the ground state, accumulated by varying the angle 
from to TT, is described by 

r d 

Pg = ~^ {il>4,\ — \ip4,)d(j), 

here is the ground state of [20| . 

Since the eigenvalues of the Hamiltonian H have not 
been affected by this unitary transformation, the eigen- 
vectors of new Hamiltonians can be obtained by act- 
ing the rotation operator on the eigenvectors of former 
Hamiltonian {H). In the other words IV'^) = gif,\ip) where 
\^) and are the eigenvectors of H and TJ^ respec- 
tively. However, the projection operator (Po (</')) of the 
new Hamiltonian {H^) and the projection operator (Pq) 
of unrotated Hamiltonian [H] are related by 

In the other side, the ground state of the renormal- 
ized chain will be related to the ground state of the 
original one by the transformation = Po{(t>)\'4>') ■ It is 
straightforward to show that the geometric phase in the 
renormalized chain is described by 



To investigate the geometric phase in systems, a new 
family of Hamiltonians have been introduced that can 
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FIG. 2: (Color online) The first derivative of geometric phase 
for different system size. For the limit of large system (high 
RG step), the nonanalyticity behavior of the first derivative 
of GP is obtained through the diverging. 



FIG. 3: (Color online) Scaling of the position {Xmux) of 
for different length chains. Xhiax goes to Ac as the size of the 
system increase as Ac = XMax + A''""'^"**. 



where /3g is geometric phase at the nth step of RG and 

7° is defined by ■ The expression for 7" is similar to 

7° where the coupling constants should be replaced by 
the renormalized ones at the corresponding RG iteration 
(n). In this approach, geometric phase at each iteration 
of RG is connected to its value after a RG iteration by Eq. 
([3]). This will be continued to reach a controllable fixed 
point where we can obtain the value of the geometric 
phase H,!!!. 



V. NUMERICAL RESULTS: SCALING 
PROPERTIES OF GEOMETRIC PHASE 

In this section we discuss the numerical results of our 
model. We present the evolution of (5g under RG steps 
versus A in Fig. ([T|). In the n RG step the expression 
given in Eq. ^ is evaluated at the renormalized coupling 
given by the n iteration of A given in Eq. ([2]). The zero 
RG step means a bare two-site model, while in the first 
RG step the effective two-site model represents a four-site 
chain. Generally, in the n RG step, a chain of 2""'"^ sites 
is represented effectively by the two sites with renormal- 
ized couplings. All curves in Fig. ([T|) have a kink at the 
critical point, Ac = 1, for large systems. Comparing Fig. 
P]) with our previous results [^-|8| and the results of exact 
solutions 0, , shows that application of QRG to man- 
ifest the GP properties, is quantitatively more accurate 
than its appUcation on quantum information properties. 
At the critical point, correlation length is infinite and 
fluctuations occur on all length scales which means that 
the system is scale- invariant. The non- analytic behaviour 
is a feature of second-order quantum phase transition. 
It is also accompanied by a scaling behaviour since the 



correlation length diverges and there is no characteristic 
length in the system at the critical point. 

Zhu has verified that the GP of ground state in the 
XY model in the transverse field obeys scaling behavior 
in the vicinity of a quantum phase transition [2lj. In 
particular he has shown that the geometric phase is non- 
analytical and its derivative with respect to the magnetic 
field diverges at the critical point. As we have stated, a 
large system, i.e. N = 2"+^, can be effectively describe 
by two sites with the renormalized coupling of in the n- 
th RG step. The first derivative of GP is analyzed as 
a function of magnetic field at different RG steps which 
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FIG. 4: (Color online) Scaling the maximum of -^j- for various 

size of system, and the maximum diverges as -j^\xMax ~ 
0.393 In iV. 
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lation length almost covers the size of the system ^ r-^ N 
which results in the following scaling form 



N 



FIG. 5: (Color online) Finite-size scaling for different lattice 
sizes through the RG treatment. The curves which correspond 
to different system sizes clearly collapse on a single curve. 



manifest the size of the system. We present in the Fig. 
^ the derivative of GP with respect to the coupling con- 
stant dPg/dX. It shows a singular behavior at the critical 
point as the size of system becomes large. This singular 
behavior is the result of the kink of the GP at X — Xc- 

It is found that the position of the maximum of d(3g/dX 
(XMax) tends towards the critical point like 



A. = A 



Ma 



-1/9 



where 6 c± 0.957 (Fig. 



Moreover, we have derived 

the scaling behavior of ^^Iam^^ versus IniV. This quan- 
tity is shown in Fig.(|3|), which behaves linearly and the 
scaling behavior is obtained as 



dX 



with K = 0.393. 

The exponent 6 is directly related to the correlation 
length exponent (v) close to the critical point. The cor- 
relation length exponent gives the behavior of correlation 
length in the vicinity of Ac, i.e., ^ ~ (A — Xc)^"^ ■ Under 
the RG transformation, Eq. ([2]), the correlation length 
scales in the nth RG step as 5^"^ ~ (A„ — Xc)~'^ = S./n'j^, 
which immediately leads to an expression for I^^Ue in 
terms of u and ns- Here ub counts the number of sites in 
each block. Dividing the last equation to ^ (A — Ac)"" 
gives I^^Uc N^^" , which implies 6 = l/v, since 

'j\\>'Max \^J\\k the critical point. We should 
note that the scaling of the position of maximum, Xmux 
(Fig. ([3]), also comes from the behavior of the correla- 
tion length near the critical point. In the large system 
size limit, by approaching to the critical point the corre- 



XMax+N-^/". 



To obtain the finite-size scaling behavior of -^Uajoxj 
we look for a scaling function in such away that all graphs 
tend to collapse on each other under RG evolution which 
results in a large system. This is also a manifestation 
of the existences of the finite size scaling for the GP. We 

have plotted l-ea:p(^-^|A„„J versus 7V(A-Amqx) 
in Fig. ([5]). The lower curves which are for large system 
sizes clearly show that all plots fall on each other. 

The similar scaling behaviours and their relation to 
correlation length exponent have been reported in our 
previous works [3, fl] . We have studied the static proper- 
ties of the ground state entanglement and low energy 
state dynamics of entanglement of ITF model by RG 
method respectively. These facts strongly imply the im- 
portant relation between quantum entanglement and ge- 
ometric phase, and provides an possible understanding of 
entanglement from the topological structure of the sys- 
tems. This point can be understand by noting that both 
of them are connected to the correlation functions, and 
are connected directly to each other by the inequality 

M. 



VI. SUMMARY 

To summarize, we have implement the idea of renor- 
malization group (RG) to study the geometric phase of 
Ising model in transverse field. In order to explore the 
critical behavior of the ITF model the evolution of geo- 
metric phase through the renormalization of the lattice 
were examined. In this respect we show that the RG 
procedure can be implemented to obtain the GP of a 
system and its finite size scaling in terms of the effective 
Hamiltonian which is described by the renormalized cou- 
pling constants. The phase transition becomes significant 
which shows a diverging behavior in the first derivative of 
the geometric phase. This divergence of GP are accom- 
panied by some scaling behavior near the critical point 
where the size of the system becomes large. The scaling 
behavior characterizes how the critical point of the model 
is touched as the system size is increased. It is also shown 
that the nonanalytic behavior of GP is originated from 
the correlation length exponent in the vicinity of the crit- 
ical point. This shows that the behavior of the GP near 
the critical point is directly connected to the quantum 
critical properties of the model. We get the properties of 
GP for a large system dealing with a small block which 
make it possible to get analytic results. 

We have also implemented the idea of the QRG to 
the anisotropic Heisenberg model XXZ in transverse 
field with Dzyaloshinskii-Moriya interactins to obtain the 
QCPs. We have been able to obtain the scaling behav- 
ior of the GP. This justifies our main idea and shows 
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that this method is universal and eould easily applied to comments, 
obtain the quantum properties of a large variety of quan- 
tum critical systems for general cases, especially for the 
models which are usually not exactly solvable. 
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